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Using a framework based on the 1 + 3 formalism we carry out a study on axially and reflection 
symmetric perfect and geodesic fluids, looking for possible models of sources radiating gravitational 
waves. Therefore, the fluid should be necessarily shearing, for otherwise the magnetic part of the 
Weyl tensor vanishes, leading to a vanishing of the super-Poynting vector. However, for the family 
of perfect, geodesic fluids considered here, it appears that all possible cases reduce to conformally 
flat, shear-free, vorticity-free, fluids, i.e Friedmann-Roberston-Walker. The super-Poynting vector 
vanishes and therefore no gravitational radiation is expected to be produced. The physical meaning 
of the obtained result is discussed. 
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I. INTRODUCTION 

In a recent paper [ij the 1 + 3 formalism [H has been 
used to develop a general framework for studying axially 
symmetric dissipative fluids. Besides some results ex¬ 
hibited in [l|, the above mentioned framework has been 
applied to the shear-free case As the result of such 
a study, it follows that all geodesic and shear-free fluids, 
are irrotational, and as consequence of this, also purely 
electric. Such a result holds for a general fluid (not nec¬ 
essarily perfect). 

Now, if we define a state of intrinsic gravitational radi¬ 
ation (at any given point), to be one in which the super- 
PoyntiM vector does not vanish for any unit timelike vec¬ 
tor [7H9|, then since the vanishing of the magnetic part 
of the Weyl tensor implies the vanishing of the super- 
Poynting vector, it is clear that when looking for gravi¬ 
tationally radiating sources (at least under the geodesic 
condition) we should consider shearing fluids. It is worth 
recalling that the tight link between the super-Poynting 
vector and the existence of a state of radiation is firmly 
supported by the relationship between the former and 
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the Bondi news function 0 (see m for a discussion on 
this point). 

From the comments above, and as a further step 
towards the understanding of gravitationally radiating 
sources, we shall consider in this work the simplest 
fluid distribution which we might believe to be compati¬ 
ble with a non-vanishing super-Poynting vector, namely: 
perfect fluid under the geodesic condition. 

However, our investigation shows that all possible 
models, sourced by a perfect fluid (which of course in¬ 
cludes the pure dust configuration as a particular sub¬ 
case) , belonging to the family of the line element consid¬ 
ered here, do not radiate gravitational waves during their 
evolution. 

We shall discuss about this result and its relationship 
with the fact that the process of radiation (including ab¬ 
sorption and/or Sommerfeld type conditions) is not a re¬ 
versible one. 

We shall heavily rely on the general framework devel¬ 
oped in [l|, keeping the same notation, and just reducing 
the general equations to the particular case considered 
here. These will be presented in an Appendix. 


II. THE PERFECT, GEODESIC, FLUID 

We shall consider axially and reflection symmetric per¬ 
fect fluid distributions (not necessarily bounded). For 



such a system, we assume that the line element may be 
written in “Weyl spherical coordinates”, as: 
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ds^ = -A^df+B^ {dr^ + r^de^)+C^d(l)^ + 2Gdedt, ( 1 ) 

where A, B, C, G are positive functions of t, r and 6. We 
number the coordinates = t,x^ = r, x^ = 0,x^ = (j). 

The specific form of the line © deserves some com¬ 
ments. Our original goal, here as well as in [ij and 0, 
has been to describe the gravitational radiation process, 
through the physical properties of its source. Such an 
endeavour, should, eventually, lead to the obtention of 
a specific source of gravitational radiation. With this 
purpose in mind, and in order to render the problem 
under consideration, analytically handable, we have im¬ 
posed the highest degree of symmetry compatible with 
the existence of gravitational radiation. For vacuum, it 
is represented by the Bondi metric 0- Here, follow¬ 
ing the framework developed in [l|, we have restricted 
the line element as much as possible, always allowing for 
the existence of gravitational radiation (at least for the 
most general matter distribution). As the result of these 
restrictions, we are able to manipulate the resulting ex¬ 
pressions, by purely analytical procedures. However, this 
is obtained at the price of dealing with a space-time, that 
is not the more general one, compatible with axial sym¬ 
metry (see the discussion below (1501) 1. 

The energy momentum tensor in the “canonical” form 
reads: 


Tad — {d + P)yaVp + Pgad^ ( 2 ) 

where as usual, /x, P and Vg denote the energy density, 
the isotropic pressure, and the four velocity, respectively. 
As in [T| we are assuming the fluid to be comoving in the 
coordinates of ©. 

The shear tensor is defined by two scalar functions 
CT/, ail, which in terms of the metric functions read (see 
eqs.(20-25) in Q): 


2a I + ail 


A “ cj ’ 


(3) 


2a II + ai = 


3 

^ 252^2 Q2 


AB^r^ 



C 

C 



(_A G_cy 


(4) 


where dots denote derivatives with respect to t. 

For the other kinematical variables (the expansion, the 
four acceleration and the vorticity) we have: 
the expansion 


AB^ 

7-2 

[2^ 

r^A^B^ + G2 



1 i 


A 

G 

A’^B'^ \ 


A 

+ G 



The four acceleration 

Ua = V^Va-p = aiKa + ttnLa, 


(5) 

( 6 ) 


with vectors K and L having components: 


GA'^B^r^ -L r<2 

iF„ = (0,H,0,0); L„ = (0,0,A---,0), 

(7) 

and where the two scalar functions (a/, a//) are defined 
by (see eq.(17) in i) 


ai = 


an = 


A' 

AB' 


A 


■' i g) 

A^s 

VA^B^r^ + G2 

A2 


^ A 


( 8 ) 

,(9) 


whereas the vorticity vector is defined through a single 
scalar H, given by (see eq.(29) in Q) 


H (10) 

2ABVA^B^r^ + G^’ 

where primes denote derivatives with respect to r. 

We shall restrict our system to the case of vanishing four- 
acceleration Oq, = 0. This condition implies that 


and 


a/=-^ = 0 => A = A(t,d), (11) 




or 


(f) ^ j = - J A^edt + G{r,9). (13) 

Given that G{t, 0, 0) = 0, from (IT^ we find that 

Am = 0 ^ A = A(t) and G = GA. (14) 

In this case, reparametrizing the time coordinate, the line 
element takes the form 

ds^ = —dt^ + B^ {dr'^ + r'^dd'^) + ‘2'G{r, 0)dtd0 + G'^dcjA, 

(15) 
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and the kinematical quantities become 
2B‘^r^ + &B C 

B2j.2 j^G2~b ^ C’ 


(16) 


n = 


G' 

2 B\/S 2 r 2 +G 2 ’ 


(17) 


52^2 + 2G^ B C 
S2r2 +G2 B G’ 


(18) 


O'/ — an 


3G2 B 
B^r2 + 


(19) 


Now, from the regularity conditions, necessary to ensure 
elementary flatness in the vicinity of the axis of symme¬ 
try, and in particular at the center (see 0 , 0 , 0 ), 
we should require that as r Ri 0 


= ^flW(t, 6 »)r", ( 20 ) 

n>l 


implying, because of (ITTI) that in the neighborhood of the 
center 


^ (6»)r". (21) 

n>3 


Thus, taking the time derivative of the above equation 
and combining with dm), dm) and dm> we obtain 


£ 11 —£i 


and 


O'/ — an 
0 -f O'/ 



£i 


Yt 


^{0 + ai){e + an) 
(26) ■ 


£ 11 —£ I 


ai — an 
0 -I- an 



£ii 


Yt 


i(0-kcr/)(0-kcr//) 


(27) 

where dSJ) has also been used. 

In the above expressions, £iji, are two of the three 
scalar functions defining the electric part of the Weyl 
tensor (see [ij for details) (the third scalar function is 
denoted £kl )■ Also, Yt, is one of the structure scalars 
obtained from the orthogonal splitting of the Riemann 
tensor which are defined in eqs.(38-50) in [I|.The others 
are denoted by Yjjj^KL,XT,i,ii,KL,Zijijnjv- 

From the above equations and it follows at once 
that as r ~ 0 


f/-f// = ^ [4 ”^(g0) 

n>4 




(28) 


and (1^ . unless 0 -I- c/ = 0 -I- an = 0. 

However, this last condition cannot be satisfied. In¬ 
deed the absence of singularities in 772 at r = 0 requires 
from dm, that ai ii Ri r as r ~ 0 which would produce 
0 Ri 0 , implying in turn because of dm 

yr(r = 0) = (// + 3P).=o = 0. (29) 


This last result in turn implies that as r approaches 0, 


ai — an 


^ \aP {t,e) - a^Pj^ {t,e) 

n>4 


( 22 ) 


Now, for the length of an orbit at t, 9 constant, to be 
27rr, close to the origin (elementary flatness), we may 
write, as r —>■ 0 , 


G^n{t,0), (23) 


Now, in order to satisfy (l29l) we have to assume, the 
equation of state 

+ 3P) = 0, (30) 

at least in the neighborhood of the center. Excluding 
this possibility on physical grounds, we have to assume 
Q + aiji^O. 

Next, from the combination of (lAlOl) . (IAI2I) and (IA13I) 
we obtain 


implying 


G'^^{t,e), G^e-nfi, (24) 

where ^{t,9) is an arbitrary function of its arguments, 
which as appears evident from the elementary flatness 
condition, cannot vanish anywhere within the fluid dis¬ 
tribution. 

Finally, observe that a combination of (fT 6 l) - (IT^ pro¬ 
duces 


G 

VB^r^ + G 2 


H 2 G' 

2H\/H2^2 _^G-2 


+ i7i 




= 0 , 


(31) 

where Hi,H 2 are the two scalar functions defining the 
magnetic part of the Weyl tensor. 

From the above equation it follows that either G = 
0 n = 0 , or 


H2n + Hi 


2aI -I- an 
3 


= 0 , 


(32) 


{ai-an){Br‘^+ &) = G‘^{QYai). (25) 

We shall next make use of the full set of equations 
deployed in [I|, written for the specific case considered 
here, they are given in the Appendix A. 


where dm-dsi) have been used. 

Since goes to zero at the center, as r, we are left 
with three possibilities to proceed further, namely: 

1. The vorticity is assumed to vanish (G = 0). 
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2. G ^ 0 and the term within the square bracket in 
HMD does not vanish at the center, meaning that 
we may write 

cr/ = ^ CT// = ^ 4/^^"- (33) 

n—0 n—0 

3. G 0 and the term within the square bracket in 
dsn does vanish at the center implying that 

CT/ = ^ an = (34) 

n>l n>l 

Next, contracting (A.7) in [ij with K and L we obtain, 
respectively 

P' = 0 ^ P = P{t,0), (35) 


GP^t + Pe = 0, (36) 

which, due to the regularity conditions on the axis of 
symmetry, implies that either G = 0 and P = Pit), or 
G ^ 0 and P = constant. 

We shall now analyze the three possible cases mentioned 
before. 

A. G = 0 


whereas from (IA2I1 - (IA4I) . we obtain, respectively 

d-+ ^0CT +Yf = 0, (41) 

Ykl = 0^ Xkl = 0, (42) 

and 

d-id^ + ^edd-Y//= 0. (43) 

Then, from (HII1) - (H51) we get 

y^=yn=Y=l^P^, (44) 

and using (42)-(45) and (47)-(50) in [I| it follows that 

Xi = Xii =X, £i= £ii = £, Y = -X = £, 

£kl = Xkl = Ykl = 0. (45) 

Two constraint equations follow from (IA6I) and (IA7I) 

C 

20' -a - 3d— = 0, (46) 

o 

20,e - d,e - 3d^ = 0, (47) 

whereas from (jASp and (IA9p we get 


From the well established link between radiation and 
vorticity (see [a and references therein), it might be 
inferred that no gravitational radiation is expected to be 
emitted in this case. 

In what follows we shall provide a formal proof of this 
result. 

As mentioned before, in this case we have P = Pit); 
then, from , (|^J) and we get 


0 = 2 


B C 


B 


B 


aj=an = a=- 


C 

G’ 


(37) 


implying: 


0 = 2d- 


3G W 


C 




B 


(38) 


Next, equation is the Raychaudhury equation for 

this case, which reads 


0 + ^0^ + + Yr = 0, 


(39) 


or, using 


^ 1~2 2^. 3C Yt ^ 

+_0^ + __ + __o. 


(40) 


H,=- 


- * ( 
2Br \ 

'G,e 

. c 

O' y 

\ (dG),, 

1 2BrC ’ 

(48) 

a 

^ W 

(1 

G'\ 

+ g) 

(dG)' 

2BC ■ 

(49) 


From the two equations above we find that if the mag¬ 
netic part of the Weyl tensor vanishes, then 


dG = 4’it), (50) 

where ipit) is an arbitrary integration function. The 
above equation implies, because of the regularity con¬ 
dition Cit, 0, 9) — 0, that d = 0. 

This last result in turn implies (as mentioned before) 
that our line element o (or (IT5]) ') is not the more gen¬ 
eral one, since in the spherically symmetric limit, it does 
not contain the Lemaitre-Tolman-Bondi metric. Also, 
the Szekeres metric [l6l - [l^ (its axially symmetric version) 
cannot be recovered from dTSl) . In other words, the shear 
of the fluid (in our models) is sourced by the magnetic 
part of the Weyl tensor. Therefore in the spherically 
symmetric limit, we recover the shear-free case. 

Another important conclusion emerges from (H51) , (1491) . 
Indeed, since B is regular at the origin (r Ri 0), whereas 
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C behaves as C ~ r, then for the magnetic Weyl tensor 
to be regular at the origin, we must demand there a ~ r 
(at least). But then, (1^ or (l4^ implies that £ ps r at 
the origin. 

Also, from (ITfl) . it follows that Q g ^ 0 (in the neigh¬ 
borhood of the center), and then from (l38l) it follows that 

^ = m, ( 51 ) 

7 


Now, combining the two equations above with (1581) we 
may write 

^ = %^C' = -f{r,e)Cg, (61) 

implying 

C = K(t)C(r,0), (62) 

and, because of dSZl 


where 1(0) is an arbitrary function of its argument. 
Next, using (|151) and (ITTl) we may write 


iJ2 = -r-f(r,0)Hi. 
Then, using (ICTl) in (1551) we get 


(63) 



II 

(52) 

^ = ^l(r,O) + eir,0), 

(64) 




with 


-(1) 

II 

(53) 

'1,91' + 77 'r^ — 1 
r(l -I- r^ 7 ^) ’ 

(65) 


.0 


implying 


Feeding back the two equations above into (l48)l and (l4^ 
we obtain 

R = T(t)P(r,d)^| = ^ = /(t), 

( 66 ) 


BrHi = - 




(54) 


On the other hand (lAllI) and (IA16I) . produce, respec¬ 
tively 


2 C" (Bryy 

~C 'WT)\ ’ 

(55) 

and 

(HiBrC^y + (H2BC^)fi = 0, (56) 

whereas the combination of the equations above with (H51) 
and (H5t produces 


H2^B+H2y-^-—^=r 




HiC'r + H 2 Cfi = 0, (57) 


or 




(58) 


Next, using (l371) in (1521) and (l53l) . it follows that 


b\ ^ BC^g C^g 
B j BC C ’ 

/ ,9 


(59) 


bV _ BC C' 
b) ~ ~ ~C' 


( 60 ) 


where T(t) and B(r,0) are arbitrary functions of their 
arguments. 

Using the above result in (15^ . it follows at once that 
Hi = H 2 = 0 , which as mentioned before implies a = 
0. This last result in turn implies that the expansion 
scalar only depends on t and taking into account that 
Yt = 47 r(/r -I- P), the Raychaudhury equation requires 
M _ 

Next, (IA14I) reads in this case 

if' + = 0, (67) 

which due to the fact that f « 0 at the origin, implies 
f = 0. 

Thus our spacetime is conformally flat, shear-free, and 
due to the fact that the fluid is perfect, it is a FRW 
spacetime, in agreement with the result obtained in Q. 

For the sake of completeness we shall sketch another 
proof of the above result in the Appendix B. 

Now, from the fact that our system is conformally flat 
it appears that it does not radiate gravitationally (ac¬ 
cording to the criterium commented in the Introduction). 

Indeed, the super-Poynting vector can be written 
(ec.(55) in Q) as 

Pa = PiKa + Pi I La, ( 68 ) 

where according to (56) in [ij, and (HSl) . we have for the 
two scalars defining the super-Poynting vector 

Pi = 2H2Y; Pii = -2HiY. (69) 

Thus, the vanishing of f. Hi and H 2 , as in our case, 
implies the vanishing of gravitational radiation. 
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B. G 7 ^ 0 and the term within the square bracket in 
(|32|) does not vanish at the center. 


C. G 7 ^ 0, and the term within the square bracket 
in (I32|) does vanish at the center. 


In this case, from the regularity of H 2 at the center we 
may write, in the neighborhood of r r; 0 , 


Then, it follows at once from (HH) and (UHl), close to 
the center, that 


H 2 = Y^ V", (70) 

n—0 




(79) 


in which case, ((32ll implies 

(71) 

n>l 

Replacing the two above expressions in a combination of 
(lAlOl) and (IA12I) . and using (1^ . (IA3I) . (1^ . (1^ we find 
that the lowest order of r in Hi and H 2 is raised as 

Hi = Y^ H2 = Y, (72) 

n>6 n>5 

Next, from (|A14I) we may write (close to the center) 

= £jj = J2£^fr^, (73) 

n>l n>l 

for otherwise there would be an inadmissible singularity 
in the r-derivative of the energy density at the origin. 

Feeding back the two expressions above in a combina¬ 
tion of (IA12I) and (|A13I) , it follows from the lowest order 
in r, that in the neighborhood of the center 

{p + P)ai KiO{r), (74) 


Next, from the lowest order of r in (IA2|) and (IA4I) it 
appears that 

= £ii = J24fr^, ( 80 ) 

n>l n>l 


and, (HID, dMD, dUD and (HZD, as r r; 0, produce (I22L 
(E51) and 


£KL = J2^Klr^- ( 81 ) 

n>5 


Excluding singularities of the scalars Hi,H 2 , at the 
origin, we may write: 

Hi H2 = J2 1 ^ 2 "^(82) 

n=0 n—0 


Then, looking for the lowest order of r in (|A10D ~ (IA13D 
and (IA16I) we obtain respectively: 


17 ( 0 ) - _fj 

Bifi — B 


(0) 7,e 


(83) 


which implies that. 


{fi + P)Ki 0{r), 


(75) 


r(o) H-(o)^27.e B^e. „(o) 






(84) 


which of course is impossible unless we assume, close to 
the center, the equation of state 

{^l + P) = 0. (76) 


2 1 ^ ^ B 


AO), 1,9 Bfi, 


(85) 


Excluding this possible situation from physical consider¬ 
ations, we have to require that 

ai = Y, (77) 

n>l 

which in turn implies 

(Tii = (78) 

n>l 

because of (l 2 ^ . 

But of course this contradicts the main assumption of 
this case about the nonvanishing of the term within the 
square bracket in (1321) . (Eq. (l33l) l. 

Thus we have to assume (j34L implying that the term 
within the square bracket in (13211 vanishes at the center 
as r. 


= 0 , ( 86 ) 

3h[°'> + H^°^g + H^°\^ + ^) = 0. (87) 

’ 7 B 

Next, from the lowest order in dMD it follows that in 
the neighborhood of the center 



implying right there 


7 =/(t)a( 6 »), B = g{t)(i{9). 


(89) 
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Then integrating (|83)) we obtain 

„(o) _ ^ 

' a( 0 )’ 


(90) 


+ ^) = 0 - 


( 101 ) 


where x(t) is an integration function. 

An equation derived from the combination of ()84|) and 
(EZD, can also be integrated to produce 


Then, proceeding exactly as we did before, using (1891) . 
we are lead to 




„(o) 

2 " a^0) ■ 


= constant, 


( 102 ) 


(91) 


which cannot be satisfied because of the reflection sym¬ 
metry, as argued before. 


and a combination of (IMl) and (|^ also produces 


Therefore we must put 


2h[°'> + ^ = 0 . 

B 

(92) 

^2 = E^ 2 "^ 1 '"’ 

n>2 

(103) 

From the equations above it follows that 




„x(t) P e 

— 2 —— = —= constant, 
yit) a/ 31/2 

(93) 

Hi = 

n>2 

(104) 

However, this last equation cannot be satisfied. In¬ 
deed, because of the reflection symmetry, we have that 
/3(0) = / 3 ( 7 r), implying that Pfi must have a change of 
sign in the interval [ 0 , tt], whereas a and P are positive 
defined. 

Thus we must put = 0. Accordingly we 

have: 

At this point we have to stop the procedure followed 
so far with eouations (IA10I)-(|A13I) and (IA16I). since now 
the lowest order in r in these equations, contains terms 
not including Hi and iJ 2 - 

Thus let us turn to equations (IA 6 I). (IA 8 |) (IA9I). 

From the lowest order in (jA 8 p and (jA9p we find, re¬ 
spectively 

to 

II 

(94) 

^( 1 ) _ ^( 1 ) _ 
aj -ujj - ^ , 

(105) 



and 


Hi = 

(95) 

U^i^ = q{t, 9)a, 

(106) 

n>l 


with 


Next, multiplvinff ()Al 8 |) bv 2 and subtractinsr from 
(IA14D. we obtain at the lowest order of r 

_ 2p(i) 

Q,0 n 5 

(107) 




(n) 

II'' 


(96) 


^>2 


n>2 


Using the expressions above, we are now looking for the 
lowest order of r in (IA10I) - (IA13I) and (|A16I) . we obtain, 
respectively 


__«-(!) 7.e 


(97) 


where p{t) and q(t, 9) are arbitrary functions. 

Now, by the same arguments based on the reflection 
symmetry exposed before, it is easily concluded that 
p{t) = 0 implying = 0. Using this result 

in the lowest order of r in (IA 6 I) we obtain = 0. We 
can now feed these results back into (IA 8 I) and (IA9l) . and 
look for the lowest order in r. We obtain then that 


E in) rt 

a}/r 

n>3 


S'-'”’- 

n>3 


^^1 = 2 ^'^! '' 


(108) 




Using this last result again in (IA 6 L the lowest order in 
r now implies = 0 , which in turn, because of 
and (IA4I) implies 


- ^), 

2 1 ^ B 


(99) 




(n) 

II 


(109) 


n>3 


n>3 




Using the results above we can now return to (lAlOll - 
(190) (|A13|) and (IA16I) , since now the lowest order in r, in these 





























equations, only contains terms with Hi and H 2 . Doing so 
we shall raise the lowest order in r of Hi and iJ 2 , until 
the moment when the lowest order in these equations 
contains terms without Hi and i? 2 - Then we can go again 
through the whole cycle above. Now, it is a simple matter 
to see that this procedure may be continued as many 
times as desired, to obtain that = H^^^ = = 

= a I = afj = = 0 for any value of n > 0, 

implying in turn that at the center, these quantities as 
well as their r-derivatives of any order vanish. 

Then, assuming that all relevant variables are of class 
C"^, i.e. that they equal their Taylor series expansion 
around the center, we can analytically continue the zero 
value at the center to the whole configuration and there¬ 
fore, we obtain a conformally flat and shear-free space- 
time (i.e. F.R.W.). 


III. CONCLUSIONS 

We have shown that all possible models compatible 
with the line element (fTSl) and a perfect fluid, are FRW, 
and accordingly non-radiating (gravitationally). This 
clearly indicates that, both, the geodesic and the non- 
dissipative conditions, are quite restrictive, when looking 
for a source of gravitational waves. 

Flaving arrived at this point, the relevant question is: 
does this result make sense from the physical point of 
view? 

To answer to such a question, let us first remember 
that, already in the seminal Bondi’s paper on gravita¬ 
tional radiation(see section 6 in[T3), it was clearly stated 
that, not only in the case of dust, but also in the absence 
of dissipation in a perfect fluid, the system is not ex¬ 
pected to radiate (gravitationally) due to the reversibility 
of the equation of state. The rationale supporting this 
conjecture is very clear: radiation is an irreversible pro¬ 
cess, this fact emerges at once if absorption is taken into 
account and/or Sommerfeld type conditions, which elim¬ 
inate inward traveling waves, are imposed. Therefore, 
it is obvious that an entropy generator factor should be 
present in the description of the source. But such a factor 
is absent in a perfect fluid, and more so in a collisionless 
dust. In other words, the irreversibility of the process of 
emission of gravitational waves, must be reflected in the 
equation of state through an entropy increasing (dissipa¬ 
tive) factor. 

In order to delve deeper into this question, let us in¬ 
voke here the tight relationship between radiation and 
vorticity mentioned before (see the beginning of Sec. II 
A). 

Now, the equation (IA5I) in the general (non-geodesic) 
case (Eq. (B.5) in [I|), reads 

n^sV^ + 1(20 + + aii)n + = 0, (IIO) 


which of course reduces to (jA5l) if the four acceleration 
vanishes. 

From (IA5I) it follows at once that if at any given time, 
the vorticity vanishes, then it vanishes at any other time 
afterwards. Thus we should not expect gravitational ra¬ 
diation from a physically meaningful system, radiating 
for a finite period of time (in a given time interval), for 
otherwise such a radiation will not be accompanied by 
the presence of vorticity. 

But, what happens for the perfect (non-dissipative, 
non-geodesic) fluid? 

In this latter case, the condition of thermal equilibrium 
(absence of dissipative flux) reads (see eq. (57) in Q) 

a^ = -h^r.^, (Ill) 

where F = InT, and T denotes the temperature. 

Feeding back (llll|) into (jllOjl we get 

+ 1(20 + (Ji + an + = 0. (II2) 

Thus, even if the fluid is not geodesic, but is non- 
dissipative, the situation is the same as in the geodesic 
case, i.e. the vanishing of vorticity at any given time 
implies its vanishing for any time in the future. 

This result is in full agreement with earlier works in¬ 
dicating that vorticity generation is sourced by entropy 
gradients [2^-[2^. At the same time we confirm, by 
invoking the radiation-vorticity link, the Bondi’s conjec¬ 
ture about the absence of radiation for non-dissipative 
systems. 

Finally, two comments are in order before conclnding: 

• In a recent work [2^, the role played by magnetic 
fields in the generation and survival of vorticity, 
has been brought out. This strongly suggest that 
the inclusion of magnetic fields in the discussion of 
gravitationally radiating sources, deserves further 
attention. 

• Geodesic fluids not belonging to the class consid¬ 
ered here (Szekeres) have also been shown not to 
produce gravitational radiation . This strength¬ 
ens further the case of the non-radiative character 
of pure dust distributions. 
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Appendix A: Summary of equations for the geodesic From B.3 

case 


Below we shall write the equations of the framework 
developed in [I|, for the geodesic and perfect fluid case. 
Then, equations B.l- B.18 in [I|, read, respectively 
From B.l 


1 


(cT/ — an) n + Ykl = 0. 


(A3) 


From B.4 


0 + T g iy'i T + O’//) + hr — (Al) 


12 2 

d// + -ajj + gQo-// — -CT/ (ct/ + CT//) +17/ = (A4) 


From B.2 


From B.5 


12 2 

CT/ + —ct| + g®^’/ — gO"// (ct/ + CT//) + 1/ = (A2) 


fl + -(20 + CT/ + CT//)fl = O. (A5) 

O 


From B.6 


+ G2 


From B.7 


From B.8 


From B.9 




/GG G,9\ 

, 1 J 

i G G ) 

^ 35 1 

\ / . 



^ (i3V+G^)' 
C 2 + G2^ 


-CT// 


G' 


' 

G 

2 (B2r2 + G2^ 

> 


= 0 . 




+ G2 


|(20 — aii) g + G (20 — CT//)' 


+o/ 


E£-9£ 

B G 


( 5 gV 



- O’// 


5,9 2G, 

“b” 


G 




= 0 . 




+ 0-1 


B,9 Gg 


a_ 

c 

-G 


GG' 


2 ^B'^r'^ + G2^ j G-v/ B^r^ + 


G2 


|(2ct/ + CT//) g 


(^ 



- O’// 


5,9 2G9 , / 2B 2G 

5 G I 5 G 


(A6) 


(A7) 


(AS) 


52 = 7 (o'/ + 2 ct//)' + CT/ 


2G' {Br){Bry 
~ ~ B^r^ + G2 

1 


+ an 


2\/B2r2 + G2 


^ 2(Br)(Br)' + 3GG' 

C 2 (B2r2 + G^) 

f29-f2 ^ + G'(| + g) 


(A9) 
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From B.IO 


From B.ll 


—Si + + P)'^i ~ SiSkl + (30 + <Jii — (Ji) H—^ (2cr// + aj) 


, ^ . (hi,0 + 

VB^r^ + G 2 V C J B 


_ 2 {Br){Bry + GG' 
C 2 (B 2 r 2 + G 2 ^ 


= 0 . 


(AlO) 


^-Skl + Skl (20 — ai — an) + — {£i — £h) + — ^ H[ + Hi 


VB^r^ + G2 


H2,e + H 2 


2G' 2(Br)(Br)'+ GG' 

G 2 (B 2 r 2 + G 2 

'b G 


2Ge ~ 

G B I B G 


= 0 . 


(All) 


From B.12 


—f// + “^(m + P)(^ii + S^Skl H—^ (30 + CT/ — (T//) + — (2(7/ + (T//) 


+ lfe + i725,) + 




C*/ \/B2r2 + G2 


_ G' 

G B I B G 


= 0 . 


From B.13 


— 7 ;{Si + Sii)' — -{Si + Sii)Q — 7 r{^^ + P){^i + <^ 11 ) —( 2 ct 7 / + cr/)-^ ( 2 ( 7 / + cr//) 


\/B2r2 + G2 


i?i.e + ^ I + iF2 


(Br)(Br)' +GG' 
B2r2 + G2 


= 0 . 


From B.14 




2G' (B^r^ + G^)' 


a/ _B2r2 + G2 


G 2(B2r2 + G2) 
Skl^b + + Skl 


+ •?// 


G' (5^2+ G2)' 
~C ~ 2 (B 2 r 2 + G 2 ) 


2 B.G,~ 2 B G 

“b^ + ^g b' + g 


-■^^2(ct/ + 2(7//) - 312-ffi = 
o oJd 


From B.15 


KL 



+Si 


+ Skl ^ 


B, 


B 



+ G 


1 

3\/(B2r2 + G2 
2G« , B, 


G 


B 


+ G 



1 / N ^ Stt 

+ --ffl(2CT/ + (7//) — 3nH2 = , = 

3 3\/B2r2 + G2 


(G/i + /i,6(). 


(A12) 


(A13) 


(A14) 


(A15) 
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From B.16 


{an - ai) - — Ih[ + Hi 


2C' {B^r^ + G^)' 


+ G2 


H2.0 + GH2 + H2 


Ba2C,0 ~ B 2C 
— + -^ + G - + — 


G 2(B2r2 + G2) 

= [87r(/r + P) — {£j + £//)] £1. 


(A16) 


From B.17 


Wb^t^ + G2 
+iS// 


(f/ + £ 11 ) g + G {2£ii + £ 1 ) + £i 


Gfi ~ / G 5^2 B 




G 


G B^r^+G^B^ 
1 


G 


B 


G - + 


G B2r2 + G2 B 
GG' 


a _ 

G 2(B2r2 + G2) 

47r 


--ffi 


--Bi(30 + a//-(T/)-B2ll = - ^- 

3 3\/B2r2 + G2 




(A17) 


From B.18 


\/B2r2 + G2 


3B 

G£kl + fifL 


1 ^ ^ /2G' GG' 

— F£,+&,) +& —+ 




GG' 


G' 


G 2(B2r2 + G2) 


^ + <5 ? + S 


G 


B G 


• 1 4:7r 

H 2 + o772(30 + ct/ — (T//) — Hi£l = 

o 60 


(A18) 


Appendix B: Conformal flatness implied by the 
vanishing of vorticity 

We shall here provide an alternative proof of the result 
exhibited in Sec.II about the consequence of assuming 
vanishing vorticity. 

First of all we observe that due to (HI, dSB and dun 
we may write at r ~ 0 


the equations (IB3I) and (1551) can be rewritten as: 


or 




(Bl) 


where L{ 6 ) must be a regular function of 9. 

Then, using (IBip in (l63)) (always at r ~ 0) we obtain 



H 2 = -L{e)Hi. 

Next, from (lAlOl) and (jAl2p we obtain 


-Hi,g-Hi 


2Gg Bg 


G B 


H'+Hoi 


(B2) 


(2G' {Br)' 


V G Br 


(B5) 

(B6) 

= 0, (B7) 


The equations (iBTl) and (IB8I) can be integrated to obtain 

00 

Hi = b ^ r" [an (t) cos{n9) + /3„ (t) sin(n0)] , (B9) 

n—1 
00 

H 2 = b ^ r" [—/3„(<) cos(n0) + Q!„(t) sin(n0)] (BIO) 


(B3) t. 


where coefficients an{t), /3n{t) are arbitrary functions of 


Then, introducing the auxiliary function b{t, r, 6) defined 

by 


Br = bG^, 


(B4) 


Observe that for n = 0 we have 


Bt Bt 

Hi = -^ao{t), H 2 = (BH) 
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However, in the neighborhood of r 0, we have C ~ r; 
implying that Hi and H 2 are singular at the origin, and 
therefore we must put ao(t),/3o(t) = 0. 

Next, feeding back (IB9I) and (IBIOI) into (IB3I) . we have 
a system of equations for each order of r”. Using (IB2I) 
we obtain 

ai=Pi= 0, (B12) 

where we have explicitly used the fact that L{0) is a 


regular function of 6 (at least on the symmetry axis). 

Repeating the calculations for the consecutive orders, 
it is a simple matter to check that 

ctn = Pn = 0, for any n, (B13) 

thereby proving that Hi = H 2 = 0. The remaining of 
the proof follows as in Sec. III. 
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